Wave function recombination instability in cold atom interferometers 
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Cold atom interferometers use guiding potentials that split the wave function of the Bose-Einstein 
condensate and then recombine it. We present theoretical analysis of the wave function recombina- 
tion instability that is due to the weak nonlinearity of the condensate. It is most pronounced when 
the accumulated phase difference between the arms of the interferometer is close to an odd multiple 
of 7T and consists in exponential amplification of the weak ground state mode by the strong first 
excited mode. The instability exists for both trapped-atom and beam interferometers. 
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PACS numbers: 03.75.-b,39.20.+q,03.75.Be 

Recent experimental demonstrations of miniature cold 
atom guides [l], ||, ||, ||, beamsplitters |j| |[ j], ||, and 
Bose-Einstein condensation (BEC) on a chip |£ | open the 
possibility of ultra-precise inertial and rotation measure- 
ments via cold atom interferometry. Cold atom inter- 
ferometers split the wave function of the condensate and 
than recombine it by using guiding potentials that change 
from a single well into two separate wells and back. Both 
the splitting and recombination are adiabatically slow to 
avoid excitation of unwanted modes. The phase of the 
wave function in each of the two potential wells evolves 
independently once the wells are separated far enough 
and is sensitive to its local environment. As a conse- 
quence, these two wave functions acquire a relative ac- 
cumulated phase difference A(j>. When the two potential 
wells converge back to a single well, this phase difference 
results in an interference that is used to extract informa- 
tion about the differences between local environments. 

In this paper we present theoretical analysis of the 
wave function recombination instability that is due to 
the weak nonlinearity of the condensate. It is most pro- 
nounced when the relative phase difference A(j> is close to 
an odd multiple of tt. The instability consists in exponen- 
tial amplification of the weak ground state (symmetric) 
mode by the strong first excited (antisymmetric) mode. 
We calculate the instability growth rate and present its 
dependence on the shape of the guiding potential and 
the nonlinearity for both trapped-atom and beam inter- 
ferometers. 

Evolution of the condensate in the interferometer will 
be described by a one or two-dimensional nonlinear 
Schrodingcr equation (NLSE) 
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The wave function ip in Eq. ([!]) is normalized to unity, V 
is the linear guiding potential, time is normalized to the 
characteristic eigenfrequency ujq of this potential, and the 
spatial coordinates are normalized to the characteristic 

1/2 

length ao = (h/uioM) , where M is the atom mass. 
Finally, N is the normalized nonlinearity parameter. 

In its one-dimensional form, x = x, Eq. ([!]) describes 
a trapped-atom interferometer [jlO| O. The condensate 



is tightly confined in two transverse dimensions and is 
in the lowest transverse mode of the trap. Equation (|l|) 
is the projection of the Gross-Pitaevskii equation onto 
this lowest transverse mode. The nonlinearity param- 
eter N contains the overlap integral between the three- 
dimensional wave function of the condensate and the low- 
est transverse eigenmode of the trap and is, in general, 
time-dependent. In our analysis we assume N — const. 
Generalizations to the time-dependent nonlinearity pa- 
rameter N are straightforward. 

In its two-dimensional form, x = (x, y), Eq. (Q) de- 
scribes a beam interferometer |l2|, |^|. Here the con- 
densate cloud propagates along the z-axis in the guiding 
potential V(x, y, z) that confines the condensate in the 
x-y plane. Equation (|l|) is written in a co-propagating 
frame that is moving with the condensate. The time 
dependence of the guiding potential V(x, y, t) is then pa- 
rameterized by the longitudinal velocity of the conden- 
sate v: V(x,y,t) — V(x,y,zo + vt). The conditions of 
applicability of Eq. (Fy) are 



{Mv/K) = k p > (Moj //!) 1/2 



and 



kj > 1, 



(2) 



(3) 



where M is the atom mass and I is the characteristic spa- 
tial scale of the guiding potential V(x,y,z) along the z 
axis. The first inequality means that the energy associ- 
ated with the longitudinal motion of the condensate is 
much larger than the characteristic energy Hojo of trans- 
verse eigenmodes of the interferometer. The second in- 
equality ensures that undesirable backward reflections of 
the condensate wave function from the guiding potential 
are exponentially small and can be neglected. Addition- 
ally, this inequality allows one to neglect the d 2 /dz 2 term 
in Eq. (§. 

We start our analysis from the one-dimensional version 
of Eq. (Q) with the simple model guiding potential of the 
form 



V(x,t)= [1 + (p(t) - x 2 /2f 



I 1/2 



(4) 



FIG. 1: Schematic view of the recombination region. 




FIG. 2: Condensate density before and after the recombina- 
tion. 



where f3(t) is the control parameter. A zero value of (3 
corresponds to a single potential well of the potential V . 
Positive values of (3 split the potential V into two wells 
separated by the distance d = 2(2/3) 1/ ' 2 . 

We assume that the condensate is initially in the low- 
est weakly nonlinear mode of the single-well potential 
V((3 — 0) with the value of the chemical potential fj, that 
is about 12% higher than the eigenenergy of the lowest 
linear eigenmode w\ w 1.3. This value of /i corresponds 
to the nonlincarity parameter N = 0.5 in Eq. ([!]). We 
further assume that the condensate wave function is split 
by increasing the control parameter [3 and acquires a rel- 
ative phase shift A<fi = it — 2 x 10~ 2 before the recombi- 
nation. The recombination stage is modelled by choosing 
the control parameter j3{t) to be of the form 



/?(*) =A]n[exp(-t/T) + l], 



(5) 



where A and T are constants. Equation (JsJ) describes 
two separate wells that are linearly converging toward 
each other at large negative values of time and merging 
at large positive values of t. The schematic view of the 
recombination region is shown in Fig. [I]. Figure ^ shows 
the condensate density before and after the recombina- 
tion. The initial wave function is the combination of the 
lowest weakly nonlinear eigenmodes of the left and right 
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FIG. 3: Modal decomposition coefficients versus time. Ai is 
the ground state mode and A2 is the first excited mode of the 
guiding potential V . 



potential wells with the relative phase shift 7r, i.e., the 
first antisymmetric weakly nonlinear eigenmode of the 
potential V. Because the nonlinearity is weak, this mode 
practically coincides with the first antisymmetric linear 
mode of the potential V. The fractional contributions of 
the higher antisymmetric linear modes are about 10~ 4 . 
Additionally, the input wave function contains a small 
amount of the first linear symmetric mode of the poten- 
tial V at the level of 10 -4 . The parameters used for this 
calculation are A = 3 and T = 90. 

The value of the parameter T — 90 ensures that the re- 
combination is adiabatic and hence the input wave func- 
tion should map onto the lowest antisymmetric mode of 
the single- well potential V(f3 = 0) preserving its odd par- 
ity. Figure ^| instead demonstrates that the parity is bro- 
ken and the wave function after the recombination is a 
superposition of even and odd modes. To quantify this 
statement we introduce modal decomposition coefficients 
A n (t) of the wave function ip(x, t) onto linear eigenmodes 
4> n (x, (3) of the potential V by the relations 

00 

il>(x,t) =J2Mt)<t>n{x,P(t)), 



n=0 



A n (t)= I dxip(x,t)<j> n (x,j3(t)). 



(6) 



The linear eigenmodes of the potential V <j>{x, (3) and 
their eigenfrequencies u> are solutions of the eigenvalue 
problem 

u4>{x, (3) = 0) + V(x, (3)<j>{x, (3) (7) 

and parametrically depend on the control parameter (3. 
Figure |3| shows modal decomposition coefficients |Ai| 2 
and 1 2 as functions of time on a logarithmic scale. 
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Solid lines correspond to the results shown in Fig. ||. 
The input wave function shown in Fig. corresponds to 
\A 2 \ 2 = 1 and \Ai\ 2 = e = 10~ 4 . Figure [| clearly demon- 
strates exponential amplification of the weak mode 1 at 
the expense of the strong mode 2 until mode 2 is de- 
pleted. The dotted curve in Fig. || corresponds to a run 
with the same parameters except the initial population 
of the first mode was chosen to be e = 1CP 8 . This curve 
shows that the instability exists only in some range of 
the control parameter (3. If this parameter is above or 
below certain values, the system is stable. Indeed, in the 
run with e = 10~ 8 the first mode stabilizes at the value 
of 10~ 2 after ( k 40 despite the fact that the mode 2 (not 
shown) remains undepleted. 

To get an insight into the nature of the instability, 
we fix the control parameter (3 at some constant value 
and replace Eq. (Q) by a set of coupled equations for the 
modal amplitudes A n (t) 
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Here 



Knklm = / dx4> n 4>k4>l4>n 



(8) 



(9) 



are the intermodal overlap integrals. 

Keeping only A\ and A 2 in Eq. (^) results in the set 
of two coupled equations 

i^A t = wivli + TV {k 111x \A 1 \ 2 A 1 
+2k 112 2\A 2 \ 2 A 1 + n 1X22 A\A\) , 



i— A 2 = uj 2 A 2 + N (n 2222 \A 2 \ 2 A 2 
+2k 1122 \A 1 \ 2 A 2 + kh 22 A 2 A 2 ) . 



(10) 



Introducing new variables j/i = 2KeAiA 2 , y 2 = 2Irm4iA2 
and j/3 = |^4i| 2 — \ A 2 \ 2 transforms Eq. (|To|) to the form 



dt 
d_ 

di 
d_ 

di 



yt = [-Au + NK\ + NK 2 y 3 ] y 2 , 
y 2 = [Alu - NK\ - NK 3 y 3 ] y 1; 



y 3 = -NK A yiy 2 , 



(11) 



where 



Auj — uj 2 — u>i, 

K\ = - («mi - K2222) , 

K 2 = - (itiui + K2222 - 2^1122) • 

K 3 = 2 ( K 11H + K 2222 - 6KH22) • 

K4 = 2k 1122 . 
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FIG. 4: Instability growth rate for a trapped-atom interfer- 
ometer versus control parameter (3 for several values of the 
nonlinearity N. 



Equations (|TT|) have two integrals of motion 



vl + vl + vl = 1, 



NK A yl + 2 [-Aw + NKi] y 3 + NK 2 y 2 



(13) 



where c an integration constant, and can be solved in 
terms of elliptic functions. In the following we will be 
interested in the limit when the amplitude of A 2 (the first 
odd eigenmode) is much larger than that of A\ (the first 
even eigenmode). This limit corresponds to j/1.2 — > 0, 
j/3 — > — 1 in Eq.([ll]) and yields solutions of the first two 
Eq. ( p] ) of the form yi,y 2 oc exp(7i) with 7 given by the 
relation 



7 2 = [-Au> + N(ki 122 - K2222) 
x [Alo + N(k 2222 - 3^1122)] • 



(14) 



(12) 



The instability corresponds to 7 2 > 0. 

Both the eigenfrcquencies ui and the overlap integrals 
K in Eq. depend on the shape of the guiding po- 
tential V , i.e., on the value of the control parameter (3. 
Figure || shows the instability growth rate 7 as a func- 
tion of the control parameter f3 for different values of 
the nonlinearity N. The curves were obtained by solving 
the eigenvalue problem Eq. ([?]) and then calculating the 
growth rate according to Eq. ( |l4| ) . 

For each curve plotted in Fig. [| there exists a cutoff 
value of the control parameter (3 below which the sys- 
tem is stable. These results explain saturation of the 
dotted curve in Fig. |3|. The maximum growth rate cor- 
responds to the values of (3 of the order of one. For 
large values of (3 the growth rate becomes exponentially 
small and scales as (3 oc (Aw) 1 / 2 . The asymptotic /3> 1 
corresponds to the frequency difference Aw being expo- 
nentially small and to all the overlap integrals k being 
exponentially close to each other. Analytical estimates 
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FIG. 5: Instability growth rate for a beam interferometer 
versus control parameter f3 x for several values of the nonlin- 
earity N. 



carried out in the framework of a model of two identi- 
cal square wells show that formally there exists no upper 
boundary on the separation between the wells that makes 
7 2 negative. From the practical point of view though the 
instability is quenched once the wells are sufficiently sep- 
arated because the growth rate becomes exponentially 
small. 

The preceding analysis has also been carried out with 
other shapes of the guiding potential V and gave sim- 
ilar results. Furthermore, the recombination instability 
also exists in the case of beam interferometers with two- 
dimensional guiding. Results of the previous modal anal- 
ysis are straightforwardly carried over to this case. The 
instability growth rate is still given by the Eq. (|l4| ) . The 
only difference is that one has to use the two-dimensional 
version of the eigenvalue problem Eq. (0) . 

Figure || shows the instability growth rate for the guid- 
ing potential 



V(x,y) = [bl(x,y) + bl(x,y)+b 2 z } 1 / 2 , (15) 



where 



xy, 
l. 



(16) 



Equation ( |l5| ) describes the guiding potential for the two- 
wire cold atom interferometer ]l3| ] in the limit when the 
characteristic size of the wave function is much smaller 
than the spatial scale of change of the magnetic field. The 
P x < part of the graph is the mirror reflection of the 



(i x > part and is not shown. Comparison of Fig. [| and 
^ shows that the instability exhibits similar qualitative 
behavior both in one- and two-dimensional cases. 

Results of the preceding analysis demonstrate fasci- 
nating dynamics of weakly nonlinear guided mater waves 
and offer some insight into the design of cold atom in- 
terferometers. First, the guiding potential should not be 
made "too adiabatic". Changing the value of T from 
T = 90 to T = 30 (T = 10) in Eq. (§) reduces the total 
amplification of the weak mode from about 10 6 to 10 2 
(10 1 ) while still preserving the adiabaticity. Secondly, 
the control parameter should be chosen in such a way as 
to minimize the time the wave function spends in the re- 
gion with the largest instability growth rate. Finally, the 
instability analysis results can be used to estimate prac- 
tical sensitivity of the interferometer in terms of errors 
of determining the phase difference A(f> for a given level 
of the nonlincarity. 
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